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TECHNICAL  MEMORANDUM  1594 

A  FLAT  WING  WITH  SHARP  EDGES  IN  A  SUPERSONIC  STREAM* 

By  A.  E.  Donov 


In  this  work  there  is  given  an  approximate  solution  of  the  problem 
of  a  two-dimensional  steady  supersonic  stream  of  ideal  gas,  neglecting 
heat  conduction,  flowing  around  a  tHn  wing  with  sharp  edges  at  small 
angles  of  attack.  (Determination  of  the  law  of  distribution  of  pressure 
along  the  wing,  lifting  force  ahd  head  resistance  of  the  wing.) 


PART  I 


The  problem  of  the  investigation  of  the  mechanical  action  of  a 
moving  gas  on  an  immovable  wing  appears  as  a  special  case  of  the  some¬ 
what  more  general  problem  of  the  investigation  of  the  mechanical  action 
of  a  moving  gas  on  an  immovable  fixed  wall  constraining  the  motion  of 
the  gas.  In  our  own  explanation  we  begin  with  the  formulation  of  this 
last  problem  in  which  we  confine  ourselves  only  to  the  consideration  of 
the  steady  two-dimensional  forces  of  ideal  gases  not  subject  to  the  action 
of  gravitational  forces.  In  the  plane  of  motion  of  the  gas  we  shall 
arrange  an  immovable  rectangular  coordinate  system  in  such  a  manner  that 
it  is  situated  as  in  figure  1.  We  introduce  three  functions  v,  p,  and  p 
of  the  independent  variables  x  and  y  defined,  respectively,  as  the  veloc¬ 
ity,  density,  and  pressure.  The  vector  functions  v  will  be  determined  by 
a  pair  of  scalar  functions  of  the  independent  variables.  For  these  func¬ 
tions  we  shall  agree  to  take  either  the  functions  vx,  vy  defined  as  the 

projections  of  the  velocity  of  the  axis  x  and  y,  respectively,  or  the 
functions  v  and  (J,  defined,  respectively,  as  the  absolute  value  of  the 
velocity  and  its  angle  with  respect  to  the  positive  direction  of  the  x-axis, 
measured  in  the  counterclockwise  sense.  In  what  follows  we  limit  ourselves 
to  the  consideration  only  of  flows  for  which  the  function  £  satisfies  the 
condition 


-  £< 


<  * 


(1) 


*Izvestiia-Akademia,  NAUK,  USSR,  1 959>  PP  605-626. 


2 


NACA  TM  1  yjk 


As  Is  well  known,  the  study  of  the  gas  motion  under  consideration 
leads  to  the  investigation  of  the  following  system  of  differential 
equations 


dx  dy 


(2) 


Here  k  is  the  adiabatic  exponent  (for  air  k  =  1.^05).  If  the  motion 
of  the  gas  is  constrained  by  an  immovable  frictionless  fixed  wall  in  the 
plane  XOY,  the  gas  will  be  adjacent  to  it  along  some  curve.  We  shall 
call  this  curve  the  "contour  K." 


Consider  the  unit  Vector  t  tangent  to  the  contour  K  directed  in 
such  a  manner  that  its  projection  on  the  x-axis  is  positive.  Denote 
by  Pk  the  angle  which  it  makes  with  the  x-axis.  Clearly  0^  may  be 

regarded  as  a  function  of  ttje  abscissa  x  of  that  point  of  the  contour  K 
associated  with  the  vector  t.  We  denote  this  function  by  Pfc(x)  and 

assume  that  it  is  continuous.  If  the  function  0k(x)  *s  prescribed  and, 
moreover,  the  coordinates  of  any  point  of  the  contour  K  are  given,  the 
form  and  position  of  the  contuur  is  completely  determined.  We  agree  to 
take  as  origin  the  left  edge  of  the  contour  K.  Then  the  equation  of 
this  contour  will  have  the  form 


* 


y 


tan  0jt(x)dx 


0 


(3) 
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We  can  write  this  equation  more  briefly  if  we  designate  its  right-hand 
side  by  yk(x) 

y  =  Yk(x)  (M 

Since  in  the  flows  under  consideration  the  direction  of  the  velocity  on 
the  contour  K  must  coincide  with  the  vector  "t^  the  condition  on  the  flow 
along  an  immovable  fixed  frictionless  wall  may  be  written  in  the  following 
fashion 


3  =  Mx) 


(5) 


at  y  =  yk(x) .  The  condition  (5)  must  be  added  to  the  system  of  equa¬ 
tions  (2)  as  a  qualifying  boundary  condition.  Much  work  has  been  dedicated 
to  the  investigation  of  solutions  of  the  system  (2)  subject  to  the  con¬ 
dition  (5) .  Of  these  we  are  interested  here  only  in  those  in  which  the 
flow  is  supersonic,  i.e.,  flows  at  every  point  of  which  the  following 
condition 


v  >  a 


(6) 


is  satisfied,  where  a  is  the  local  speed  of  sound 


The  Investigations  contained  in  these  works  divide  in  two  fundamental 
directions.  The  first  direction  is  represented  in  works  in  which  solu¬ 
tions  of  the  problem  are  achieved  with  the  help  of  numerical  or  graphical 
processes  permitting  the  step-by-step  calculation  of  a  system  of  parti¬ 
cular  values  of  the  desired  functions.  (Works  of  Busemann,  Kibelia, 
and  Frankl.)  The  fundamental  achievements  of  the  methods  represented  by 
these  works  consist  of  the  fact  that  by  their  use  many  actual  practical 
problems  may  be  solved  quantitatively  of  which  the  solution  by  other 
methods  would  present  great  difficulties.  In  particular  these  methods 
solve  thoroughly  corner-nonpotential  problems.  The  chief  defect  of 
these  methods  is  that  the  solutions  obtained  are  numerical  so  that  it 
is  impossible  to  obtain  a  general  qualitative  estimate  of  the  phenomena 
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under  investigation.  The  second  direction  is  represented  by  the  works 
of  Meyer,  Ackeret,  Prandtl,  and  Busemann,  which  are  confined  to  a  culti¬ 
vation  of  an  exact  theory  of  irrotational  flows.  The  results  are  based 
on  the  fact  that  in  the  case  where  vorticity  is  absent  the  character¬ 
istic  system  of  differential  equations  (2)  admit  of  integrable  combina¬ 
tions.  This  theory  leads  to  series  of  approximate  results  of  any  desired 
accuracy,  giving  a  complete  qualitative  and  quantitative  picture  of  the 
flow.  Since  our  investigation  is  mostly  connected  with  the  theory  of 
irrotational  flows  we  give  below  a  brief  introduction  to  the  fundamental 
methods  and  results  of  this  theory. 

We  introduce  the  stream  function  1/  defined  by  the  following  relations 


a* 

dx 


dy 


> 


-pvx 


(8) 


As  is  well  known  from  equations  (2),  (7),  and  (8)  the  following  relations 
follow  without  difficulty 


(9) 


y2 

T 


+ 


(10) 


Where  0,  tQ  denote  quantities  which  display  the  flow  once  and  for  all 

as  a  function  only  of  ijf.  With  the  help  of  equations  (2),  (9),  and  (10) 
it  is  easy  to  obtain  the  two  equations 


dx  dy 


(11) 
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(a2  -  vx2)^*  +  (a2  -  vy2)^  -  vxvyi  — + 
dx 


8y 


'dv. 

—1 

idx 


8y 


=  0 


(12) 


where  Q  denotes  a  quantity  defined  as 


n  =  p 


dtr 


d  In  0 


d\|»  k(k  -  1)  d^ 


(13) 


Equations  (11)  and  (12)  represent  linear  relationships  between  the  first 
partial  derivatives  of  the  functions  vx,  vy  with  respect  to  x  and  y. 

Since  every  flow  under  consideration  is  supersonic,  the  entire  region 
of  the  flow  may  be  covered  by  a  pair  of  families  of  characteristics.  The 
differential  equations  of  these  characteristics  are  obtained  easily  by 
the  use  of  equations  ( 11)  and  (12).  For  one  family  of  characteristics, 
which  we  shall  agree  to  call  the  first  family,  we  obtain  the  equations 


dy  a  mjdx 


(a2  -  v^cos^Jm-v  +  v^sin  0  cos  0  /,s\ 

d(v  cos  0)  +  m2d(v  sin  0)  =  ft  - dx 

v2cos20  -  a2 

and  for  the  other,  which  we  shall  agree  to  call  the  second  family,  we 
have  the  equations 


dy  =  mgdx 


(16) 


/  v  ,  v  (a2  -  v2cos20)m{)  +  v2sin  0  cos  0  . 

L(v  cos  0)  +  mxd(v  sin  0)  =  fl  - - — - - - -  dx  (17) 

V2cos20  -  a2 


Here  m^,  ioq  denote  the  following  expressions 


mi 


-v2sin  0  cos  0  +  a^v2  -  a2 
a2  -  v2cos20 


(18) 
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-v2sin  fl  cos  6  -  a\/v2  _  a2 

ra2  = - - - 7 - — -  (19) 

a2  _  y^cos^p 

We  now  consider  a  supersonic  stream  with  constant  hydrodynamical 
elements  (i.e.,  functions  v,  0,  p,  .  p,  a).  We  shall  call  this  flow  the 
undisturbed  flow.  The  values  of  the  functions  v,  p,  p,  a  in  the  undis¬ 
turbed  stream  will  be  denoted  by  w,  pQ,  pQ,  respectively,  and 

the  ratio  w/ao  by  M.  Since  the  stream  under  consideration  is  super¬ 
sonic,  M  >  1.  We  shall  choose  the  direction  of  the  velocity  of  the 
undisturbed  stream  to  correspond  to  the  direction  of  the  x-axis. 

We  assume  that  the  undisturbed  stream  strikes  an  immovable,  fixed, 
frictionless  wall  (contour  K) ,  inclined  in  such  a  manner  that  in  flowing 
around  this  wall  the  stream  never  detaches  from  it  and  remains  super¬ 
sonic  everywhere.  We  may  distinguish  two  cases  of  flows  of  this  type. 

Case  I.-  The  contour  K  is  situated  in  such  a  manner  that  the 
condition 


3k(0)  £  0 


(20) 


is  fulfilled.  In  this  case,  as  is  well  known,  there  appears  a  curve  of 
weak  discontinuity  OC  (figs.  2(a),  and  2(b))  proceeding  from  the  origin 
and  dividing  the  entire  flow  in  two  parts.  On  one  side  of  the  curve  of 
weak  discontinuity  OC  extends  the  region  containing  the  undisturbed 
stream  and  on  the  other  the  region  of  flow  around  the  wall.  In  the  region 
of  flow  around  the  fixed  frictionless  wall  the  hydrodynamical  elements  of 
the  stream,  generally  speaking,  are  not  constant  but  vary.  In  what  fol¬ 
lows  we  shall  call  this  part  of  the  stream  the  disturbed  stream.  In  the 
entire  region  of  the  flow  under  consideration  the  functions  v,  p,  p,  p,  a 
are  continuous  but  their  partial  derivatives  with  respect  to  x  and  y 
(all  or  only  some)  exhibit  jump  discontinuities,  at  least  on  the  curve  of 
weak  discontinuity  OC.  The  same  curve  OC  appears  as  a  characteristic  of 
the  second  family  since  the  hydrodynamical  elements  of  the  stream  are  con¬ 
stant.  On  this  line  the  following  relationships  will  hold  in  the  entire 
region  containing  the  stream 


(21) 
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9  =  0O 


(22) 


where  0q  denotes  a  quantity  defined  as 


0O  = 


P0_ 

P0k 


(23) 


From  equations  (13),  (21),  and  (22)  we  easily  obtain 


SI  =  0 


(24) 


i.e.,  the  flow  under  consideration  is  irrotational.  By  virtue  of  rela¬ 
tion  (24)  the  right-hand  side  of  equations  (15)  and  (17)  vanish  and  these 
equations  can  be  integrated.  As  a  result  of  integration  of  equation  (15) 
we  obtain  the  relationship 


p  +  <p(v)  =  constant  (25) 


satisfied  along  any  characteristic  of  the  first  family,  and  as  a  result 
of  integrating  equation  (17)  we  have 


3  -  <p(v)  =  constant  (26) 


satisfied  along  a  characteristic  of  the  second  family.  cp(v)  denotes  a 
function  defined  as 


<P(v) 


k  +  1 


arc  tan  1 


k  +  1 


-  arc  tan 


(27) 


Since  on  the  curve  of  weak  discontinuity  OC  the  quantities  v  and  p  have 
the  values  w  and  0,  respectively,  the  following  relation  is  satisfied 
along  every  characteristic  of  the  first  family  intersecting  this  line 
and  consequently  in  the  entire  region  of  the  disturbed  stream: 
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P  +  q>(v)  =  cp(w) 


(28) 


From  equations  l6,  26,  and  28  it  immediately  follows  that  the  char¬ 
acteristics  of  the  second  family  (the  curve  OC  being  among  these)  are 
straight  lines  since  along  each  of  these  characteristics  the  hydrodynam- 
ical  elements  are  constant. 

Making  use  of  these  circumstances  it  is  not  difficult  with  the  aid 
of  equations  28,  26,  22,  21,  l6,  10,  7>  and  5  to  construct  expressions 
for  the  functions  v,  p,  p,  p,  a  in  the  region  of  the  disturbed  flow. 
However,  the  construction  of  these  expressions  is  not  of  great  interest 
since  our  chief  interest  is  centered  on  the  construction  of  an  expres¬ 
sion  for  the  pressure  on  the  contour  K  which  may  be  accomplished  without 
the  use  of  these  expressions  for  the  hydrodynamical  elements  of  the  flow. 
Actually  equation  28  allows  us  to  determine  the  velocity  v  as  a  function 
of  the  angle  of  inclination  of  this  velocity  with  the  x-axis  at  every 
point  of  the  region  filled  by  the  disturbed  flow.  Since  by  virtue  of 
equation  5  the  angle  of  inclination  of  the  velocity  with  respect  to  the 
x-axis  is  a  given  function  of  x  on  the  contour  K  there  is  the  possibility 
of  using  equations  22,  21,  10,  9 1  and  7  to  determine  the  pressure  p  as  a 
function  of  x  on  the  flow  around  a  contour.  If  we  limit  ourselves  to 
the  consideration  of  slightly  disturbed  flows,  i.e.,  flows  whose  hydro- 
dynamical  elements  differ  but  little  from  the  hydrodynamical  elements  of 
the  undisturbed  flow,  the  expression  for  the  pressure  on  the  flow  around 
a  contour  K  may  be  written  in  the  form  of  a  series.  This  series  has  the 
form 


P  =  Po  +  q  aipk(x)  +  a2Pk2(x)  +  a^p^x)  +  a1+pk4(x)  +  .  . 


(29) 
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where 


q  ■ 


PQw2  P0kM2 


ax  =  2(M2  -  l)'1/2 


a2 


.  (M2  -  1)  ”^f2  -  3i2  +  M4^ 


aj  =  (M2  -  1) 


-7/2 


—  -  2M2  +  ^(k  +  1)M4  +  ~  ^  +—  M6  +  ^-±-i  M 


8 


a*  -  (M2  -  1)  -  f  M2  ♦  M4  +  + 


25. 


+  30k  -  Sk^  -f  3k^  m8  +  -21  -  20k  +  3k2  +  2k?  M10  +  ?  +  gfc.rJg  M12) 
12  48  48  / 


Case  II.-  The  flow  around  a  contour  K  is  situated  in  such  a  manner 
that  the  following  inequality  is  satisfied 


pk(0)  >  0  (30) 


In  this  case,  as  is  well  known,  a  line  of  shock  discontinuity  OD  appears 

(fig.  3)  proceeding  from  the  origin  0  and  dividing  the  entire  flow  under  j 

consideration  in  two  parts.  On  one  side  of  this  line  is  the  region  of 

the  undisturbed  stream  and  on  the  other  the  region  in  which  the  fluid 

flows  around  the  fixed  frictionless  wall.  Just  as  in  case  I  we  call  the 

flow  in  the  region  in  which  the  stream  around  the  fixed  frictionless  wall 

is  accomplished  the  disturbed  flow.  In  the  present  case,  in  contrast  to 

case  I,  the  functions  v,  p,  p,  p,  a  exhibit  Jump  discontinuities  on  the 

shock- line  OD. 


10 


NACA  TM  13^4 


In  the  region  of  the  disturbed  flow  these  functions  must  satisfy 
not  only  equations  2,  5,  and  7  but  also  the  dynamical  conditions  across 
the  shock  line.  Considering  the  flow  to  be  only  slightly  disturbed,  these 
conditions  may  be  written  in  the  following  form 


,  a 


2  w2  .  a0£ 


2  k-1  2  k  -  1 


(31) 


v  =  v(l  +  b]_p  +  b2P2  +  b^3  +  bi^  +  .  .  .) 


(32) 


where 


bx  =  -(M2  - 


1 

1)'  2 


b2  =  -(M2  -  l)-2(i  +  M4 


b*  a  -(M2  1)"  2  7  +  -  M2  +  2(k  -  1)M^  +  #  ■  ~  ^  1  t-1).  m8 

6  2  4  24  32 


b4  *  -(M2  -  1)“5 


(i.  +  |„2  +  dLt 

\24  8  24 


2Tk  ±  12k2  ^ 
24 


3  +  3k  -  k2  +  K?  m5  +  -3  -  k  -  3k2  +  m3jq\ 
16  48  / 


_0_ 

e0 


1  +  +  .  .  . 


(33) 


r  * 
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where 


i3  »  Ife-i.:..!!  m6(m2  -  i)"  i 
J  12 


i4  .  m^(m2 

12 


d4 


+  2(k  -  2)M2  -  (k  - 


dx 


eo  +  ex0  +  e^2  + 


(54) 


where 


1 

e0  =  (M2  -  1}"  2 


el  = 


H-if^M2  -  if2 
4 


Condition  31  shows  that,  disregarding  the  presence  of  Jump  discon¬ 
tinuities  in  the  functions  v,  {3,  p,  p,  a,  equation  21,  Just  as  in  case  I, 
is  valid  throughout  the  entire  region  filled  by  the  flow  under  considera¬ 
tion.  However,  condition  22  is  not,  generally  speaking,  fulfilled  in 
the  case  now  under  consideration.  However,  there  is  the  possibility  of 
speaking  of  satisfying  this  condition  approximately.  In  fact,  consider 
equation  33.  Its  right-hand  side  does  not  contain  terms  in  the  first 
and  second  powers  of  0.  Therefore,  for  slightly  disturbed  flows,  equa¬ 
tion  22  may  be  regarded  as  approximately  satisfied  on  the  line  OD  and  con¬ 
sequently  throughout  the  entire  region  filled  by  the  flow  under  considera¬ 
tion.  From  this  it  follows  that  in  the  region  of  disturbed  flow  equa¬ 
tion  24  may  be  regarded  as  approximately  satisfied,  which  means  that 


4 
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* 

where 


l*  =  -ft2  ~  1)  M^(M2  -  1)"  2 
?  12 

u  =  M^(M2  -  1)"5U  +  2(k  -  2)M2  -  (k  -  1) 

12 


dy 

dx 


=  eo  +  e-^p  +  e2P2 


+  •  •  • 


(54) 


where 


1 

e0  =  (M2  -  1)"  2 


-  l)'2 
4 


Condition  31  shows  that,  disregarding  the  presence  of  jump  discon¬ 
tinuities  in  the  functions  v,  p,  p,  p,  a,  equation  21,  just  as  in  case  I, 
is  valid  throughout  the  entire  region  filled  by  the  flow  under  considera¬ 
tion.  However,  condition  22  is  not,  generally  speaking,  fulfilled  in 
the  case  now  under  consideration.  However,  there  is  the  possibility  of 
speaking  of  satisfying  this  condition  approximately.  In  fact,  consider 
equation  33.  Its  right-hand  side  does  not  contain  terms  in  the  first 
and  second  powers  of  p.  Therefore,  for  slightly  disturbed  flows,  equa¬ 
tion  22  may  be  regarded  as  approximately  satisfied  on  the  line  OD  and  con¬ 
sequently  throughout  the  entire  region  filled  by  the  flow  under  considera¬ 
tion.  From  this  it  follows  that  in  the  region  of  disturbed  flow  equa¬ 
tion  24  may  be  regarded  as  approximately  satisfied,  which  means  that 
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equations  25  and  26  hold  on  characteristics.  For  values  of  p  and  v 
near  0  and  w,  respectively,  equation  28  may  be  written  in  the  form 
of  a  series 


v  =  v(l  +  bx'p  +  b2'p2  +  bj'p5  +  b^'p^  +  .  .  .)  (35) 


where 


1 

b-L '  =  -(M2  -  1)’  2  =  b± 


b2'  =.  -(M2  -  l)'^i  --  b2 


bj'  =  -(M2  -  1)’  2  \  +  -  M2  +  2(k  -  ljM4  +  ~~~ m6 

6  2  4  12 


b4*  ,  -(m2  -  1  )*[±-  +  —  M2  +  ^  +  LzJP*.  t  l6k2  + 

*  \24  8  24  24 


3  -  2k  -  5k2  +  4k3  „a  .  -3  +  8k  -  7k2  +  2k^  wlD\ 

r  ,  ,  M  t  ^  M  I 

32  9 6  / 


Comparing  equations  32  and  35  ve  see  that  for  slightly  disturbed 
streams  the  first  may  be  substituted  for  the  second  with  good  approxi¬ 
mation.  Consequently,  for  slightly  disturbed  flows,  equation  28  will  be 
approximately  satisfied  along  the  line  OD.  Since,  on  the  other  hand, 
along  each  characteristic  of  the  first  family  equation  25  is  approximately 
satisfied,  equation  28  will  be  approximately  satisfied  throughout  the 
entire  region  of  disturbed  flow.  The  approximate  expressions  for  the 
functions  v,  p,  p,  p,  a  are  constituted  exactly  like  the  accurate  expres¬ 
sions  for  these  functions  in  case  I.  Substituting  the  approximate  expres¬ 
sion  for  the  function  p  in  the  right-hand  side  of  equation  $k,  we  obtain 
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a  differential  equation  of  the  first  degree  for  the  approximate  deter¬ 
mination  of  the  form  of  the  shock  line.  Summing  up  our  considerations 
we  can  deduce  that  the  accurate  results  contained  in  case  I  can  serve  as 
approximate  results  for  case  II,  and  further  that  expression  29  can  serve 
as  an  approximate  expression  for  the  pressure  on  the  flow  around  a  contour 
in  case  II.  These  same  considerations  show  that  there  is  no  sense  in  cal¬ 
culating  all  terms  in  this  expression.  It  is  sufficient  to  limit  ourselves 
to  the  first  two  or  three  terms. 

From  all  that  has  been  said  about  cases  I  and  II  one  may  conclude 
that  the  form  of  the  contour  K  may  be  made  up  in  such  a  manner  that  art¬ 
fully  constructed  shocks  may  be  caused  to  appear  in  the  region  of  flow 
around  the  fixed  frictionless  wall.  In  such  cases  when  we  pay  attention 
to  this  phenomenon,  the  results  we  have  obtained  are  valid,  not  for  the 
entire  region  of  flow  around  the  fixed  frictionless  wall,  but  only  for 
that  part  in  the  neighborhood  of  the  front  side  of  the  flow  around  a 
contour.  The  fundamental  problem  of  the  present  work  is  the  construction 
of  approximate  expressions  for  the  pressure  on  the  flow  around  a  contour 
in  case  II,  with  the  calculation  of  the  circulation  of  the  flow  occasioned 
by  the  presence  of  the  shock  discontinuity  OD.  In  spite  of  the  fact  that 
in  the  case  of  the  presence  of  circulation  it  is  impossible  to  integrate 
equations  15  and  17,  there  is  the  possibility,  however,  of  making  up  such 
combinations  of  differentials  from  equations  14,  15,  16,  and  17,  adding 
to  these  equations  expressions  for  differentials  of  the  stream  function, 
that  with  the  aid  of  these  combinations  it  is  possible  to  construct  expres¬ 
sions  which  we  shall  integrate.  Investigations  concerning  the  preceding 
construction  constitute  the  contents  of  the  following  section. 


PART  II 


Suppose  we  have  a  flow  corresponding  to  case  II  of  the  preceeding 
section.  Assume  that  in  this  flow  the  hydrodynamical  elements  in  the 
region  of  the  disturbed  stream  differ  infinitely  little  from  the  hydro - 
dynamical  elements  in  the  region  of  the  undisturbed  flow.  We  revamp 
somewhat  our  notion  of  the  region  of  disturbed  flow.  Shortly  before  we 
agreed  to  apply  this  name  to  the  region  bounded  by  the  curvelinear 
triangle  made  up  of  the  curve  OC2  (contour  K),  the  shock  line  OC]_,  and 

the  characteristic  of  the  first  family  C1C2  emerging  from  the  lowest 

point  of  the  contour  K  (fig.  4).  Taking  into  consideration  equations  5, 
14,  18,  and  34,  it  Is  not  difficult  to  conclude  that  with  the  assumptions 
made  just  now  relative  to  the  hydrodynamical  elements  the  curvelinear 
triangle  OC1C2  differs  infinitely  little  from  the  isoceles  straight-line 

triangle  O'Ci'Cp'  (fig-  5)  vhere  the  equal  sides  O'C^'  and  Ci'Cp'  are 
parallel  to  characteristics  of  the  second  and  first  families  in  the 
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undisturbed  flow.  As  for  the  functions  ^(x),  0,  v,  p,  p,  a  we  assume 

that  they  all  have  the  properties  of  differentiability  and  continuity  to 
as  many  degrees  as  may  be  necessary  to  insure  legitimacy  of  operations 
which  are  performed  upon  them.  Moreover,  we  assume  that  in  the  flow 
under  consideration  the  infinitesimal  quantities  0^(x),  0jt'(x),  p^'Cx), 

0,  v  -  w,  p  -  p q,  p  -  pQ,  a  -  a 0  have  the  same  order  of  magnitude. 

Taking  this  last  group  of  infinitesimals  as  fundamental  (having  unit 
order  of  magnitude)  we  shall  agree  in  what  follows  to  adhere  to  the  fol¬ 
lowing  system  of  notations  appearing  in  investigations  involving  infinitely 
small  quantities.  By  £m  (m  being  any  positive  integer)  let  ua  denote 

an  infinitesimal  whose  order  of  magnitude  is  not  less  than  m.  Clearly 
such  a  mode  of  notation  does  not  exclude  the  possibility  of  several  dif¬ 
ferent  infinitesimals  being  denoted  by  the  same  symbol,  and,  vice  versa. 

The  same  infinitesimal  may  be  denoted  by  several  different  symbols.  Thus, 
for  example,  if  an  infinitesimal  a  is  denoted  by  e^,  the  infinites¬ 
imal  2a  may  also  be  denoted  by  e^,  and,  moreover,  the  infinitesi¬ 
mals  a  and  2a  may  be  denoted  by  €j,  ei* 

On  an  arbitrary  characteristic  of  the  second  or  first  family  the 
equation 


dijr  -  pv(sin  0dx  -  cos  0dy) 


(56) 


will  be  satisfied  by  virtue  of  equation  8  throughout  the  entire  region 
filled  by  the  flow.  Eliminating  dx  and  dy  from  equations  14,  15, 
and  36  and  taking  into  account  formulas  13  and  21,  we  arrive  at  the 
equation 


d(v  cos  0)  +  mgdCv  sin  0)  *  Gjd  In  0  (37) 

which  is  satisfied  on  any  characteristic  of  the  first  family.  Here 
denotes  the  quantity 

a2  v^sin  0  cos  0  -  (v2cos20  -  a2)m-j| 

-  - - - - -  (58) 

k(k  -  I)v(v2cos20  -  a2)(m^  cos  0  -  sin  0) 
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On  the  other  hand,  having  the  integral  25  of  the  equation 

d(v  cos  p)  +  m2d(v  sin  p)  =  0  (39) 

it  is  easy  to  find  an  integrating  factor  Lq  of  this  equation,  such  that 
after  multiplying  by  Lq  it  may  be  written  in  the  form 

d(p  +  *(v)]  =  0  (40) 

In  order  to  determine  L-^  we  have  the  obvious  relationship 

LqJtL(v  cos  p)  +  m2d(v  sin  p)J  =  d|p  +  cp(v)J  (41) 

from  which  we  obtain  without  difficulty 

Ll(®2  cos  P  ”  sin  P ) vcLJ3  =  dp  (4-2) 


i  consequently 

) 

I 

Lq  =  - - -  (43) 

*  v(m2  cos  p  -  sin  p) 


If  now  we  multiply  both  sides  of  equation  37  by  Lq,  this  equation  takes 
the  form 

d[p  +  9(v)]  =  Hqd  In  9  (44)  I 

i 

| 

I 

where  Hq  denotes  the  quantity 


(v2cos2p 
Hq  =  - 


v^sin  p  cos  p 


(45) 


k(k  -  l)v2 
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We  denote  by  Hiq  the  value  of  at  v  =  w,  p  =  0.  We  have 


H10  -  ,  -  -  -  1)1/2 

xv  k(k  -  1)M2 


(46) 


Equation  44  may  be  rearranged  In  the  following  fashion 


d[p  +  <p(v)l  =  H10d  In  i-  +  (Hi  -  Hio)d  In  i-  (47) 

0O  90 


Now  choose  an  arbitrary  point  S  in  the  region  of  disturbed  flow  and  lead 
a  characteristic  of  the  first  family  through  it.  We  denote  the  point  of 
intersection  of  this  characteristic  with  the  shock  line  by  A  (fig.  6). 
Integrating  both  sides  of  equation  47  along  the  above  characteristic  from 
point  A  to  point  S  we  obtain 


0s  +  <P(vs)  -  0a 


where  06,  v8,  0B  denote,  respectively,  the  values  of  0,  v,  0  at 
the  point  S  and  0a,  va,  0a  denote  the  values  of  these  quantities  at 
the  point  A.  Taking  account  of  equation  (52)  we  have 


va  =  v(l  +  b1pa  +  bj@a2  +  bjSa3  +  b4Pa**  +  .  .  .) 

(49) 

We  introduce  the  quantity  vn-|  defined  with  the  help  of  the  expansion  55 
in  the  following  fashion 

Va!  -  w(l  +  bjj}a  +  b20a2  +  bj’pii5  +  V0a^  +  •  •  •)  (cf  eq.  (55)  -Tr.) 

(50) 
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By  this  definition  of  the  quantity  we  have 

Pa  +  <P(val)  =  9(v)  (51) 

With  the  help  of  formulas  49,  50,  and.  51  we  rearrange  the  expression 
Pa  +  <P(va)  111  the  following  manner 

Pa  +  <P(va)  =  Pa  +  <P(vai)  +  <P(va)  -  <P(vai) 


*  <p(v)  +  <p(va)  -  <p(val) 

»  9(w)  +  <p'(w)[(va  -  w)  -  (VgX  -  w)J  + 

i<p"(v)[(va  -  w)2  -  (vjq  -  w)2J  +  .  .  . 

-  <p(w)  +  wp'(w)£(bj  -  bJ,)Pa3  +  (ty.  -  H')pa^j  + 
v^>”(w)b1(b3  -  bj'JPa1*  +  e5  (52) 


Calculating  cp'(w),  <p"(w)  we  obtain 


<p'(w)  - 

wbi 

(53) 

2bo 

(54) 

<p"(v)  -  — 

w2b^3 
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Using  formulas  52,  53,  and.  54  we  easily  find 


Pa^  +  e5 


(55) 


Now  pass  a  stream  line  through  the  point  S  and  denote  by  P  the  point  of 
intersection  of  this  line  with  the  shock  line.  Since  the  stream  function  C 
is  constant  along  this  line  we  have 


3a  +  <P(va)  =  q>(w)  -  — (b5  -  b5')pa5  + 


2b  p 


(b5  -  b3')  -  — (b4  -  b4') 


08  = 


(56) 


where  0p  denotes  the  value  of  0  at  the  point  p.  Taking  logarithms  of 
both  sides  of  equation  (33)  we  obtain 

In  i-  -  1^3  +  l4'p4  +  i5«p5  +  .  .  .  (57) 


Since  the  values  of  the  coefficients  l4',  ,  .  .  .  will  not  be  needed 

in  what  follows,  we  shall  not  calculate  them. 

Using  formulas  (56)  and  (57)  we  easily  see  that 


In  -  In  =  l5(Pp5  -  pa3)  +  I4'(pp4  -  pa4)  +  e5  (58) 

eo  eo 


where  0p  denotes  the  value  of  0  at  the  point  p.  Assuming  that  the 

mean  value  theorem  is  applicable  to  the  integral  arising  from  the  right 
hand  side  of  equation  (48),  we  easily  find! 


■^Instead,  take  a  slightly  more  general  assumption  admitting  the  part 
AS  of  the  characteristic  under  consideration  to  be  divided  in  the  same 
finite  number  of  parts  in  such  manner  that  on  each  part  the  mean  value 
theorem  can  be  applied  to  the  integral  under  investigation. 
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/  (Ha.  -  Hio)d  In  i-  =  l3(H!  -  Hio)(3p3  -  0&5)  +  €  (59) 

J  AS  e0  ^ 

(in  consequence  of  this  equation  one  must  keep  in  mind  that  frj.  -  H^q  =  €^) . 
Here  denotes  the  value  of  at  some  point  on  the  characteristic 

under  consideration  between  the  points  A  and  S. 

Using  relations  (55)  >  (58)  >  and  (59)  we  write  equation  (48)  in  the 
following  form 


0s  +  <P(vs)  =  <p(w) 


A(bj  -  b3')pa5  + 
bl 


Mbk  -  b4') 
bl 


+ 


Hio^(0p3  -  0a3)  +  Hlol4'(0pJ+  -  Pa4)  + 

l5(Hi  -  Hlo)(0p3  -  0a3)  +  e5  (60) 


We  denote  by  B  the  intersection  of  the  characteristic  of  the  first 
family  under  consideration  with  the  contour  K.  Applying  formula  (60)  to 
the  point  B  (which  is  possible,  since  the  point  S  was  chosen  arbitrarily) 
we  obtain 


0b  +  9(vb)  =  $(*)  _  j~(b3  '  b5' )Pa3  + 


2b, 


bl 


2<b5  .  b5')  -  ~(b4  -  b4') 


0a4  + 


Hlob^O3  -  Pa3)  +  Hi0l4'(P04  -  Pa4)  + 
lj(Hi  -  H1O)(0O3  -  Pa3)  +  e5 


(61) 
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where  pb,  v^  denote,  respectively,  the  values  of  p  and  v  at  the 
point  B  and  Po  denotes  the  value  of  (3  at  the  point  0. 

We  now  proceed  to  the  derivation  of  an  expression  for  pa.  From 
formula  (6o)  we  have 


0s  +  <P(vs)  ■  <J>(v)  +  €j  (62) 

From  equation  (62),  using  formula  (35)  we  obtain 

v8  -  w(l  +  b!0s  +  b2082)  +  ^3  (63) 

and  denoting  by  mgg  the  value  of  nig  at  the  point  S  we  obtain,  by 
using  formulas  (19)  and  (63) 

m2s  -  (M2  -  l)"1/2  +  il-iM^M2  -  l)"2  ps  +  e2 

*  e0  +  2eiPs  +  £2  (&0 

Analagous  to  the  derivation  of  equation  (47) ,  which  holds  on  character¬ 
istics  of  the  first  family,  we  may  derive  equation 


djp  -  q>(v)J  ■  HgQd  In  ~  +  (H2  -  HgoJd  In 


_9_ 

©o 


(65) 


which  is  valid  on  characteristics  of  the  second  family.  Here  H2  denotes 
the  function  defined  as 


H2 


(v^cos^p  -  a2)m2  -  v^sin  p  cos  p 
k(k  -  ljv2 


(66) 


and  H20  denotes  the  value  of  this  function  at  p  =  0, 


v  =  w. 
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Now  pass  a  characteristic  of  the  second  family  through  the  point  S 
and  denote  by  Q  its  intersection  with  the  contour  K.  Integrating  both 
sides  of  equation  (65)  along  this  characteristic  from  the  point  Q  to  the 
point  S  we  obtain 


-  <P(vq) 


(67) 


where  vq,  0q  denote  respective  the  values  of  p,  v,  0  at  the 

point  Q.  Since  the  contour  K  is  a  stream  line  we  have 

0q  =  0<°)  (68) 


where  6^°)  denotes  the  value  of  0  at  the  point  0.  Assuming  that  the 
mean  value  theorem  can  be  applied  to  the  integral  arising  from  the  right- 
hand  side  of  equation  (67)  we  easily  find,  with  the  aid  of  formulas  (56), 
(57) »  and  (68) 


Ps  “  Pq  “ 


<p(v8) 


€3 


(69) 


Applying  formula  (62)  at  the  point  Q  we  have 

Pq  +  <p(vq)  -  <P(v)  +  ej  (70) 


Eliminating  cp(w)  from  equations  (62)  and  (70)  we  arrive  at  the  fol¬ 
lowing  equation 


Ps 


<P(vq) 


e3 


(7D 


22 


naca  m  1394 


Families  (69)  and  (71)  give 


Ps  =  Pq  +  e3  (72) 

On  the  shock  line  we  take  an  arbitrary  point  F  (fig.  7)  and  pass 
through  it  a  characteristic  of  the  second  family  in  the  region  of  the 
disturbed  flow  and  we  denote  by  pf,  m 2f,  respectively,  the  values  p 

and  m2  at  the  point  F.  Applying  formula  (64)  at  the  point  F  we  obtain 

m2f  =  e0  +  +  £3  (73) 

We  denote  by  (  —  ;  the  slope  of  the  tangent  to  the  shock  line  at  the 
VcLx/f 

point  F.  Fron  equations  (34)  we  have 


— i  =  eQ  +  e^Pf  +  €2  (7^0 

ax  f 


Comparing  formulas  (73)  and  (74)  we  see  that  the  characteristic  of  the 
second  family  passing  through  F  and  the  shock  line  at  this  intersection 
make  an  infinitesimal  angle  with  each  other  moreover,  if 


Pf  >  0 


(75) 


the  slope  of  the  characteristic  of  the  second  family  is  greater  than  the 
slope  of  the  shock  line  at  the  point  F.2 


2It  is  easy  to  show  that  if  the  shock  line  is  unbroken  and  moreover 
condition  (30)  is  satisfied  the  inequality  p  <  0  is  impossible  on  this 
line.  As  a  matter  cf  fact,  in  the  opposite  case  the  shock  line  is  broken 
since  with  p  <  0  condition  (34)  must  be  replaced  by  the  following  con¬ 
dition  in  virtue  of  Tsemplen's  theorem 


dy 

dx 


"  (e0 


e^P  + 


0 
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Denoting  by  L  the  intersection  of  the  characteristic  of  the  second 
family  under  consideration  with  the  contour  K  and  by  x^  abscissa  of  this 

point,  we  have 


*Z  =  *1  (76) 

Let  Pz  denote  the  value  of  p  at  the  point  L.  Using  equations  (5) 
and  (76)  and  MacLauren's  formula  we  obtain 

Pi  -  Pk(°)  +  Pk'(°)xi  +  €3  (77) 

Applying  formula  (72)  at  the  point  F  we  obtain 

Pf  =  Pi  +  €3  (78) 

As  a  consequence  of  equations  (77)  and  (78) 

pf  =  pk(0)  +  pk* (0)xI  +  e3  (79) 


Since  the  point  F  was  chosen  arbitrarily  on  the  shock  line  by  use  of 
equations  (74),  (76) ,  and  (79)  we  can  obtain  the  following  differential 
equation  for  the  shock  line 


—  =  e0  +  ejpjJO)  +  €2 


(80) 


Consequently  the  equation  of  the  shock  line  may  be  written  in  the  fol¬ 
lowing  form 


y  = 


+  eiPk(°) 


x  +  e2 


(81) 


24 


NACA  TM  1394 


Applying  formulas  (64)  and  (72)  to  an  arbitrary  point  situated  on  the 
characteristic  of  the  second  family  FL  we  easily  obtain  the  differential 
equation  of  this  line  from  the  following  form 


—  =  e0  +  2e-]J32  +  e2  (82) 

dx 

Employing  formulas  (76)  and  (77)  this  equation  may  be  written 

—  53  eQ  +  ^iP^O)  +  €2  (83) 

dx 

Consequently  the  equation  of  the  characteristic  FL  may  be  written  in  the 
form 

y  -  yi  +  [e0  +  26^(0 )J(x  -  XX)  +  e2  (84) 

where  y^  denotes  the  ordinate  of  the  point  L. 

On  the  other  hand,  taking  account  of  formulas  (3)  and  (76)  we  have 


Employing  formulas  (85)  and  (76) 


tan  pk(x)dx  -  €2  (85) 

we  may  write  equation  (84)  in  the  form 


-  gqX-i  +  e2  (86) 

Applying  formulas  (8l)  and  (86)  at  the  point  F  and  denoting  by  Xf,  yf 
the  coordinates  of  this  point  we  obtain 


4J 
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* 


yf  =  [e0  +  elPk(°) 


Xf  +  e. 


yf  =  |  eQ  +  2eipk(0) jxf  -  e0x2  +  e2 
From  equation  (&[)  we  easily  obtain 


C1 

XZ  =  —  XfPk(O)  +  e2 
e0 


(87) 


(88) 


Replacing  Xj  in  the  right-hand  side  of  equation  (79)  by  the  expression 
in  formula  (88)  we  obtain 


Pf  -  Pk(°)  +  XfPk^'Pk*  +  e3 
e0 


(89) 


We  denote  by  xa,  ya  the  coordinates  of  the  point  A  and  by  xb,  y^ 

the  coordinates  of  the  point  B.  Applying  formula  (89)  at  the  point  A 
we  arrive  at  the  following  result 

Pa  =  Pk(°)  +  ~  xaPk(0)Pk’(°)  +  e3  (90) 

e0 


We  now  express  xa  in  terns  of  x^.  To  this  end,  using  formulas  (l4) 

and  (18),  we  write  the  differential  equation  for  characteristics  of  the 
second  family  in  the  following  fashion 


dy 

-  =  -e0  + 


(91) 
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Employing  formula  (91)  we  write  the  equation  for  the  characteristic  AB 
of  the  first  family  in  the  form 

y  =  Yb  -  eo(x  -  Xb)  +  €X  (92) 

Taking  account  of  formula  (3)  we  have 


Yb 


Pxb 

/  tan  Pk(x)dx  = 

J  0 


consequently  equation  (92)  may  be  written 

y  =  -e0(x  -  xb)  +  ex 


(93) 


m 


On  the  other  hand,  equation  (8l)  for  the  shock  line  may  be  written  in 
the  form 


y  =  e0  x  +  ex 

and  applying  formulas  (94)  and  (95)  at  the  point  A  we  obtain 

\ 

Ya  ■  -eo(*a  ”  xb)  +  el 

> 

Ya  -  eoxa  +  «i 

Fran  equation  (96)  we  easily  find 


(95) 


(96) 


x 


a  ~ 


+  €1 


(97) 
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And  consequently 


Pa  =  Pk(°)  +  xb3k(O)pk'(0)  +  e. 

2e0 


(98) 


Substituting  this  expression  for  pa  in  the  right-hand  side  of  equa¬ 
tion  (6l)  and  substituting  Pq  for  Pk(°)  111  the  fundamental  formula  (5) 

we  obtain 


Pb  +  ^(vb)  *  <P(v)  -  — (b3  -  b3')pk3(o) 


l  2b  p 

— (b4  -  b4') - -(bj  -  bj') 


bl 


Pk  (°)  - 


*1 

[b3 .  bj- 

2e0 

1 — 1 
& 

_ 1 

Hiob 


xbPk5(°)Pk'(0)  +  e5  (99) 


Employing  relation  (35)  >  ve  easily  obtain  from  equation  (99) 

/ 

vb  =  w  *  1  +  bxpb  +  b2Pb2  +  bj’pb5  +  b4'pb4  +  (b3  -  b3')pk3(0)  + 


2b- 


b4  -  b4' - (b3  -  b y) 


Pk^(0)  +  ^(b3  -  b3')pk5(o)pb  + 


3e 

2ei 


i[b3  -  bj'  +  HloIjb1JxbUlt3(O)0k,(O) 


+  ec 


(100) 
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Substituting  v^,  (3^,  x^  for  v,  £k(x)  and  x,  respectively,  in 

fomula  (100)  we  arrive  at  the  following  final  expression  for  the  velocity 
on  the  contour  K: 

v  =  w  •  1  +  bxBk(x)  +  b2pk2(x)  +  b5'pk5(x)  +  bi/p^U)  +  (bj  -  b5')pk5(0)  + 

2bp  1  2b o  - 

b4  -  b4'  — ~(b3  “  b3'^  (°)  +  “j^(b3  "  b3'  )Pk^(°)0k(x)  + 

3^  r 

~j[b3  -  b3'  +  HloZ3byx0k5(o)pk'(O)  ■  +  e5  (101) 
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On  the  other  hand,  by  virtue  of  equations  (5)  and  (33)  the  following 
equation  holds 


s(0)  ,  1 

- -  1  +  1^(0)  +  *4Pk4(0)  +  .  .  .  (105) 

e0 


Substituting  the  expressions  for  v  and  obtained  in  formulas  (101) 

and  (105),  respectively,  in  the  right-hand  side  of  equation  (104)  we  obtain 
after  elementary  transformations  the  desired  formula  for  the  calculation 
of  the  pressure  on  the  contour  K 


p  =  p0  +  q  aipk(x)  +  a2Pk2(x)  +  aj0k5(x)  +  aj^U)  +  aldpk3(0)  + 


a2dPk+(°)  +  a3dPk5(°)Mx)  +  a4dPk5^0)pk'(°)x] 


+  €r 


(106) 


where 


aid  =  -2(b5  -  bj')  - 


21 5 


k(k  -  1)M2 


=  LIImV  - 1,-  f(.  i  + 


»2d  58  _2(b4  -  b4*)  +  — (b3  -  b3*) 


21, 


k(k  -  1)M£ 


m4(M2  -  i)-5/l  *+*  +  **  ?k  -  M2  +  :1Q  -  a  ±&LzJL  ^  + 
1  2  4  8 


9  -  Vi2  +  2k3  ^  t  -3  +  k  4-  3k2  -  k^  mq’ 
16  32 
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,  0  4bp\  2l3bx 

a3dL  =  (b3  "  b5* )  (2M2bl  -2b! - -  )  + 


bl/  k  -  1 


t^(M2  -  1)_5|-  +  7,  +  ^  ~  5k  m2  + 

6  24 


•4  +  3k  +  6k2  -  k3  ^  +  3  -  7k  -  7k2  +  yk?  ^ 


24 


96 


a4d  =  — ~(b3  ■  b3'  +  Hl0*3bl) 


=  (k*-.1)2  M8(m2  -  if  5(-l  + M2  +  J\ 

16  \  2  8  / 


For  x  =  0  the  formula  (106)  takes  the  form 

P  =  P0  +  q[aiPk(0)  +  a2pk2(0)  +  a5'pk3(0)  +  a4,pk4(0>]  +  €5  (107) 


where  af  =  a^  +  ald,  a4'  =  a^  +  a2d  +  a^. 

Formula  (107)  may  be  used  for  the  calculation  of  the  pressure  on  a 
flat  plate  which  is  inclined  at  an  angle  Pk(0)  to  the  undisturbed  flow. 

In  order  to  single  out  of  the  right-hand  side  of  equation  (106)  those 
terms  which  depend  exclusively  on  the  presence  of  the  shock  in  front  of 
the  contour  K,  we  add  to  the  contour  K  under  consideration  an  arc  O'O 
of  finite  length  in  such  a  manner  that  this  arc  is  tangent  to  K  at  the 
point  0  and  is  parallel  to  the  x-axls  at  O'  (fig*  8).  Since  the  flow 
around  such  an  additional  contour  is  accomplished  without  the  appearance 
of  shocks  (we  suppose  that  the  angle  between  the  direction  of  flow  and  the 
x-axis  and  the  derivative  of  this  angle  with  respect  to  x  are  both 
infinitely  small),  formula  (29)  may  be  employed  in  the  calculation  of  the 
pressure  on  this  contour.  Comparing  formulas  (29)  and  (106)  and  denoting 
by  £pstoss  Pressure  resulting  from  the  presence  of  the  shock  front, 

we  obtain 
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^stoss 


aldPk"^®)  +  a2dPk*(®)  +  a3dPk^^^k^x^  + 


HdPk^0)^0^ 


+  e5 


(108) 


We  may,  in  turn  single  out  of  the  expression  for  A?stoss  the  "t®1111 

depending  solely  on  the  vorticity  caused  by  the  presence  of  the  shock. 

In  order  to  do  this  we  add  to  the  contour  K(0C2)  a  straight-line  segment 

tangent  at  the  point  0  (segment  O'O  in  fig.  9)*  With  the  contour  0'0C2 
a  shock  is  formed,  but  the  shock  line  O'C^  is  straight  so  that  vortex 
formation  is  absent.  Calculating  the  pressure  on  the  portion  OC2  of  the 
contour  0 '0C2,  we  obtain 

p  -  Pq  +  qfaiPjtCx)  +  a2Pic2(x)  +  a^pk^(x)  +  a^p^x)  + 


aldPk3(°)  +  a2dPk^°)  +  a3dPk5(°)Mx^]  +  e5  (109) 


Comparing  formulas  (106)  and  (109)  and  denoting  by  tipTot  the  pressure 
due  to  vortex  formation  earned  by  the  shock,  we  obtain 

4?rot  =  9a4dPk^0)Pk'(°)x  +  €5  (n°) 


PART  III 


We  now  apply  the  results  obtained  to  the  calculation  of  the  lifting 
force  and  head  resistance  of  a  flat  wing  with  sharp  front  and  rear  edges 
placed  in  a  supersonic  stream  having  constant  hydrodynamical  elements. 

We  place  the  origin  0  at  the  front  edge  of  the  wing  and  arrange  the 
coordinate  system  so  that  the  positive  x-axis  corresponds  to  the  direction 
of  the  velocity  of  the  undisturbed  flow  and  measure  angles  in  the  manner 
used  heretofore.  Segment  OC2  connecting  the  front  and  rear  edges  (fig.  10) 
will  be  called  the  chord  of  the  wing  as  in  the  theory  of  wihgs.  The  length 
of  this  curve  will  be  denoted  by  T  and  the  angle  it  makes  with  the  x-axis 
by  p. 
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The  form  of  the  wings  we  are  investigating  is  defined  by  a  pair  of 
contours  like  that  investigated  in  the  preceding  section,  possessing  a 
pair  of  common  points  0,  C2»  Comparing  ordinates  of  points  on  these 
contours  having  the  same  abscissa,  we  call  the  upper  contour  Ku  that 

contour  of  which  every  point  on  the  ordinate  is  greater  than  the  corre¬ 
sponding  point  on  the  ordinate  of  the  other  contour,  which  we  call  the 
lower  contour  K}.  The  function  p^(x)  for  the  upper  contour  we  denote 
by  pku(x)  and  for  the  lower  by  Pfci(x). 

We  choose  an  arbitrary  point  A  on  the  chord  of  the  wing  and  denote 
the  distance  OA  by  t.  Through  A  we  pass  a  straight  line  perpendicular 
to  the  chord  of  the  wing  and  denote  by  Ay  and  A i,  respectively,  the 

intersections  of  this  straight  line  with  the  upper  and  lower  contours. 

With  the  point  Au  we  associate  a  unit  tangent  vector  ty  and  at  the  point  Aj 

a  unit  tangent  vector  tj.  The  vectors  ty  and  t i  will  be  directed  in  such 

a  manner  that  their  projections  on  the  direction  OCg  are  positive.  We 
denote  by  pu  and  Pj,  respectively,  the  angles  these  vectors  make  with 
the  vector  OC2.  Clearly  pu,  Pj  may  be  regarded  as  functions  of  t. 

We  denote  by  pu0  and  PiQ  the  values  of  pu  and  Pj  at  the  point  0 
and  the  values  of  the  derivatives  of  pu  and  Pj  with  respect  to  t 
at  the  point  0  by  PyQ1  and  P2o'‘ 

The  abscissa  x  of  Au  may  be  calculated  from  the  formula 


_  r  t 

x  =  t  cos  p  -  sin  p  /  tan  pudt 

J  0 


OH) 


and  that  of  the  point  Aj  from  the  formula 


x  *  t  cos  p  -  sin  p  /  tan  Pjdt 


L 


(112) 


The  value  of  the  functions  pku(x)  at  the  point  Ay  is  determined  by  the 
relation 


Pku^  =  P  Pu 


(113) 
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and  the  value  of  the  function  pjc^(x)  at  the  point  Aj  by  the  relation 


PklM  =  P  +  Pi 


(114) 


Moreover  we  have  the  relations 


tan  pu4t  =  0 


♦ 

(H5) 


tan  0jdt  =  0 


(U6) 


We  assume  that  0  and  also  0U  and  0^  and  other  derivatives  with 

respect  to  t  are  infinitesimal,  quantities.  From  equations  (115)  and 
(116)  we  easily  obtain 


i  ♦  e5 

I  /0v« + b 


(117) 

(118) 


Proceeding  now  to  the  calculation  of  the  lifting  force  and  head  resistance 
of  the  wing  under  consideration,  we  remark  that  on  the  top  side  of  the 
wing  a  shock  appears  when  and  only  when 

P  +  Puo  >  0  (H9) 


and  at  the  bottom  side  when  and  only  when 


P  +  Pio  <  o 


(120) 
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We  introduce  the  quantities  a^u,  a£u>  a^u,  a^u  defined  as  follows 

alu  =  ald 
a2u  =  a2d 

if  P  +  PuQ  >  ® 

a3u  =  a5d 
a4u  =  a4d 

alu  -  a2u  »  a3u  =  a4u  =  0  lf  P  +  ho  -  0  t122) 

In  an  analogous  way  we  define  the  quantities  a-^,  &2.i>  &$l>  a]^ 

all  a  ald 
a21  -  a2d 

if  P  +  Pl0  <  0 

a3l  ■  a3d 
a4Z  "  a4d 

all  3  a2l  3  a3l  3  a4l  3  ®  if  P  +  Pio  ^  0  (i24) 

Denoting  by  Py  the  pressure  on  the  upper  contour  and  by  pj 
the  pressure  on  the  lower  contour  we  easily  obtain,  with  the  help  of 
formulas  (106),  (121),  (122),  (123),  (124) 

Pu  3  PO  +  4  [flPlcuM  +  a2Pku2^x^  +  a3Pku5M  +  a4Pku4(x)  + 

aluPku3(°)  +  a2uPku4(°)  +  aJuPku3(0)eku(x)  ♦  a4uPta3(°>»ku’  (°>x]  +  e5 

(125) 
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n  =  PO  +  ^[a^Cx)  +  a2(3kl2(x)  -  ajP^x)  +  a4Pkl^(x)  - 

allfta5(°)  +  a2i&ki^°)  +  a3iPkl5(°)Pkl(x)  +  a4iPki^0^kl‘(0)x]  +  e5 

(126) 


Let  P  denote  the  resultant  vector  of  the  hydrodynamlcal  force 
acting  on  a  unit  length  of  the  wing  under  consideration.  We  have  then 


P  = 


(12?) 


where  n  denotes  a  unit  vector  normal  to  the  contour  of  the  wing  and 
directed  inwards. 


We  introduce  the  dimensionless  coefficient  of  the  lifting  force  Cy 
and  the  dimensionless  coefficient  of  head  resistance  Cx.  These  coeffi¬ 
cients  are  defined  by  the  formulas 


(128) 


(129) 


-4 

where  Py  and  Px  denote  the  projections  of  the  vector  P  on  the 
x-  and  y- axe s,  respectively.  From  formulas  (127),  (128),  (129)  we  have 


Cy  =  -  —  f1  Pudx  +  i  [\i*x 
y  qT  J  0  qT  J  0 


(130) 


36 


NACA  TM  1394 


1  f  ^  pi 

°X  =  +  qT  J0  PU  ^  PkuU)dX  "  J0  Vl  t8n  Pki(x)djc 

where  l  denotes  the  abscissa  of  the  point  C2-  With  the  aid  of  formu¬ 
las  (111),  (112),  (113),  (114),  (117),  (118),  (125),  (126),  (130), 
and  (131)  after  a  few  elementary  transformations  we  obtain 

Cy  =  Cyl  +  Cy2  +  CyJ  +  Cylj,  +  €5  ( 132) 

where 

cyl  3  "SaiP 

°y£  “  -  j  f*  (pu2  -  Pi2)« 

Cy3  =  (ax  -  2a5)p5  -  alu(p  +  p^)5  -  an(g  +  0lo)3  + 

i(ai  -  3a3)p  (pu2  +  Pz^idt  +  -  ajj  f  (pu5  +  Pl5)dt 

4  4  -  »  » 

Cy4  =  "a2u^P  +  PuO)  +  a2l($  +  Pio)  “  a3uP(P  +  PuO)^  + 

a3lP(P  +  Pio>3  -  |  a4u(P  +  Puo)3Pvc'  +  |  a4l(0  +  Plo)5Plo'  + 


-  pz2)dt  +  a2  -  4a^p 


(Pu5  -  Pz5)dt 
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Ox  =  CX2  +  Cx3  +  Cx4  + 


(133) 


where 


-2  al 

CX2  =  2a]p  +  — 


-T  /  \ 

(pu2  +  Pl2Jdt 


cx3  a 


rQT  (pu2  -  Pi2)dt  +  y  Iq  (Pu3  '  Pl3)dt 


cx4  *  fsa 3  -  P  +  aiuP(P  +  P*>)’  +  ailP(P  +  Pio)^  + 


f  P  -  t)p'2  /0  (p»2  +  p*  jdt  + 


|^3  -  yjp  fQ  p  +  Pz3jdt  + 


?(aj + 7)  /o' 


Pu4  +  Pi  )d^ 


Then 


Let  us  consider  a  numerical  example. 

Suppose 

k  =  1.405,  M  ■  1.5,  Pq  3  1*033  kg/cm2 

q  *  —  =  1.633  kg/cm 2 

ax  =  2(M2  -  l)-1/2  =  I.789 
a2  =  (M2  -  1)“2(2  -  342  +  1.203M4)  »  2.296 
aj  =  (M2  -  1)"T/2(i.333  _  2M2  +  4.00^  -  I.815M6  + 
0.400&8)  »  3*082 

a^  »  (M2  -  1) “^(0.6667  -  O.6667M2  +  5.6l6M^  -  3.824m6  + 
2.96SM8  -  0.7840M10  +  0.0799a*12)  =  8.290 


(134) 


>  (135) 


(Equations  continued  on  next  page). 
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ald  =  1.203M4(M2  -  l)'7^ -0.3333  +  0.265^2  -  0.0327^)  =  0.2766  ] 

a  2d  =  M^(M2  -  1)'5(-1.203  +  1.317M2  -  0.643:^  +  Q.Qk’i'&fi  + 
O.0U855M8)  *  0.MJ48 

ajd  »  f^(M2  -  ir5(-0.4008  -  0.002^2  +  0.3863^  -  0.128^) 
a^d  =  0.36l5M8(M2  -  1)“5(-1  +  0 . 797^l2  -  0.09813^)  =  0.9035 

Let  us  take  as  the  functions  pu,  pj 


>(1X  ) 
f  leone  .1 


=  0.3318 


Pu  =  -2p  +  M 

T 

Pi  -  0 


Moreover  ve  assume  that 


p  <  0 


(136) 


(137) 


The  form  and  position  of  the  profile  of  the  wing,  determined  by  equa¬ 
tions  (136)  and  condition  (137) »  is  shown  in  figure  11.  It  is  easily 
3een  that  the  straight  line  drawn  perpendicular  to  the  wing  through 

its  mid  point  is  the  axis  of  symmetry  of  the  profile  under  consideration. 
Prom  equations  (136)  and  condition  (137)  we  have 


P  +  PuO  =  -P  >  0  I 


P  +  Pio  a  P  <  0 


(138) 
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Consequently, 


alu  =  aLZ  =  aid  =  0-27^6 


a2u  =  a2 l  =  a2d  ■  0.4448 


a3u  ■  a3Z  0  a3d  =  0-3318 


(139) 


a4u  =  a4l  ■  a4d  -  0.9035  I 


Using  formulas  (132),  (133),  (135)>  (136),  and  (139)  we  obtain 


Cy  a  -2aqP  -  j  a202  +  -  6a p^  + 


r10  56  _  .  \-4 

a2  -  —  a^  +  2a^d  +  2a1+djp  +  €5 


-3-578p  -  3-o6lp2  -  14.32P5  -  82.73^^  +  .  .  . 


Cx  -  y  alP2  +  4a2p5  +  / ^  +  ^  a^jp4  +  €5 


■  5-963P2  +  9-l84p3  +  4o.8p^  +  .  .  . 


I  (140) 


Let 


. 

T  =  100  cm 


(141) 
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Then 

Cy  =  0.2936 

► 

Cx  =  0.04168 

Knowing  Cy,  Cx,  T,  and  q,  we  easily  obtain 

Py  =  qTCy  =  47.9  kg /cm 

► 

Px  85  qTCx  =6.81  kg/cm 

swim 


(142) 


(143) 


In  the  present  work  the  problem,  of  a  flow  of  stream  of  ideal  gas 
around  a  thin  wing  at  small  angles  of  attack  is  investigated,  this  stream 
being  supposed  to  be  two-dimensional,  stationary,  supersonic  and  deprived 
of  heat-conduction. 

In  the  initial  part  of  the  work,  the  problem  is  stated,  and  the  well- 
known  results  obtained  by  Ackeret,  Prandtl,  and  Busemann  are  cited.  These 
results,  as  known,  are  obtained  on  the  basis  of  the  potential  supersonic 
streams  theory,  which  is  founded  on  the  existence  of  integrable  combina¬ 
tions  of  characteristics  of  differential  equations  concerning  this  problem, 
and  in  which  some  peculiarities  of  the  dynamical  conditions  on  the  line  of 
the  shock  wave  are  utilized. 

In  the  second  part  the  approximate  solution  of  the  problem  is  given 
with  an  allowance  for  vortex-formation  caused  by  the  change  of  entropy 
along  the  shock  wave,  when  receding  from  the  leading  edge  of  the  wing, 
near  which  this  shock  wave  is  formed.  For  this  purpose  differential 
equations  of  characteristics  non  admitting  integrable  combinations  are 
to  be  dealt  with.  The  solution  is  obtained  by  means  of  a  special  method, 
which  enables  us  to  find  the  approximate  integrable  combinations  of  dif¬ 
ferential  equations  of  the  characteristics.  The  obtained  combinations 
let  us  receive  the  approximate  formula  of  pressure  in  any  point  of  the 
contour  of  the  wing  investigated.  From  this  formula  the  term  is  easily 
segregated  depending  exclusively  on  the  vortex  formulation,  caused  by 
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the  change  of  entropy  along  the  shock  wave.  The  characteristic  dis¬ 
tinction  of  this  term  of  the  obtained  formula  of  pressure  from  the  other 
ones,  is  that  it  includes  the  curvature  of  the  wing  contour  at  the  leading 
edge  and  the  distance  from  this  edge  up  to  the  element  of  the  wing  for 
which  the  pressure  is  calculated. 

In  the  third  part  of  the  work  the  expressions  for  lift  and  drag 
coefficients  of  the  wing  are  given,  on  the  base  of  the  formula  of  pres¬ 
sure  obtained  above.  In  conclusion  a  numerical  example  is  studied. 


Translated  by  R.  Shaw 
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Figure  2(b). 
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Figure  6. 


